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ABSTRACT 

The energy dependence of the anisotropy (the anisotropy energy spectrum) of the large-scale diffuse 
gamma-ray background can reveal the presence of multiple source populations. Annihilating dark 
matter in the substructure of the Milky Way halo could give rise to a modulation in the anisotropy 
energy spectrum of the diffuse gamma-ray emission measured by Fermi, enabling the detection of a 
dark matter signal. We determine the detectability of a dark-matter-induced modulation for scenarios 
in which unresolved blazars are the primary contributor to the measured emission above ^ 1 GeV and 
find that in some scenarios pair-annihilation cross sections of order the value expected for thermal 
relic dark matter can produce a detectable feature. We anticipate that the sensitivity of this technique 
to specific dark matter models could be improved by tailored likelihood analysis methods. 
Subject headings: dark matter, galaxies: active, Galaxy: structure, gamma rays: diffuse background 
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1. INTRODUCTION 

The matter density of the universe is dominated by 
non-luminous, non-baryonic dark matter, though its na- 
ture remains one o f the greatest outstanding questions in 
physics (see, e.g.. iBergstroirJ [20001 : iBertone et afl 120051 
for reviews). Many weak-scale candidate dark matter 
particles can pair-annihilate to produce gamma-ray pho- 
tons and other standard model particles, providing a 
means of indirectly detecting dark matter and probing 
its prop erties. The Fermi G amma-ray Space Telescope 
[Fermi] lAtwood et al.l [20091 ). which is currently moni- 
toring the gamma-ray sky from 20 MeV to more than 
300 GeV, has generated intense interest in the possible 
detection of a dark matter annihilation signal. 

Primary targets for indirect dark matter searches 
in gamma rays include the M ilky W ay halo (e.g . , 
iSpringel et al|[2008l: iBaltz et al.li 2008: Fi eri et al.ll2009t). 
its substructure (e.g., iPieri et al.. ■2008t ~ lKuhlen et ahl 



[2008t [BaETeFaLl [200I iPieri et al.|[2q0£), the Galactic 



Cent e r (e.g., iBerezinskv et al.' '19 94|: Dodelson. et al.l 
[2OOI iSerpico fc Zaharijas 2008; Bahz et al.l l2008D. and 
nearb y Milky Way sate lli te ga la xies (e.g.. IBaltz et al.l 
2OOOI: [Tasitsiomi et al.l [200l iStrigari et al.l 120081: 



ulations of structure formation predict that the dark 
matter halo of the Milky Way is po pulated by an 
enormous number of den se subhalos (|Diemand et al.l 
[2008: Sm- ingel et al.l [20081 ) . and theoretical calculations 
indicate that cold dark matter could form structures 
on scales of ~ or smaller, far be low those that can 
currently be resolve d in simulations ([Green et al.l 120051 : 
iProfumo et"aI1 120061 : iBringmannl 120091 ). Outside of the 
Galactic center region, most photons from Galactic 
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Baltz et all 120081 : lAbdo et al.l r2010b[ ). Numerical sim- ^ 



dark matter annihilation are expected to originate in 
substructure. Almost all subhalos are likely to be too 
faint to be resolved individually by Fermi, and instead 
will contribute to the measured diffuse emission. The 
radial number density distribution of the substructure 
is not very centrally concentrated, since subhalos near 
the Galactic center are easily destroyed by tidal effects. 
Consequently, on large angular scales, the gamma- 
ray flux frorn sub s tructure appears almost isotropic 
(IKuh len et all[200l [Siegal-Gaskinsl[2?M ISpringel etal] 
2008). On smaller angular scales, the dumpiness of 
the Galactic dark matter d i stribu ti on would i nduce 
anisotropics (lSieg:al-Gaskinsl [200l iLee et al.l 120091: 
iFornasa et all 120091: 1 Andoll2009[ ). 

We explore the possibility of detecting an annihilation 
signature from the substructure of the Milky Way dark 
matter halo in the large-scale isotropic gamma-ray back- 
ground (IGRB). In this work, we use the IGRB to refer 
to diffuse emission of extragalactic and possibly Galactic 
origin which appears isotropic on large angular scales, 
but may contain fluctuations on small angular scales. 
The IGRB is thought to be composed largely of gamma 
rays from unresolved blazars (active galaxies with jets 
aligned with the line of sig ht: iSt eckcr fc Salamon' 199 _ 
Narumoto fc Totanil 120061 : iDermer ,200 7; Vent ers et a! 
2009f ) and photons produced from the interaction be- 
tween cosmic rays and diffuse gas in normal galax- 
ies ([Pa vlidou fc Fields 2002; Thomps on et al.l 120071 : 
lAndo fc Pavlidoul 120091 : IFields et al.l l20To| ) with possi- 
ble contributions from other sources such as Galactic 
and e x tragalactic dark matte r annih i lation (lUllio et al.l 



2005 



200l lElsasser fc Mannheim [2005t lOda et al . 

Zavala et al.l l2009t lAbazaiian et al.l 120101 : lAbdo et al.l 
2010aD 



To extract information about dark matter from the 
IGRB, it is necessary to disentangle the dark matter 
annihilation signal from gamma rays arising from other 
sources. One way to do this is to examine the anisotropy 
of the gamma-ray flux on small angular scales (~ few 
degrees). Anisotropics in the IGRB are also expected 
in the collective emission from unresolved members 
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of extragala ctic source classes (lAndo fc Komatsu 
Ando et al.ll2007l: iMiniati et al.ll2007t iCuoco et al 



2006 



2008 



Taoso et a ljl2009l: iFornasa et al.ll2009l: lAndo fc Pavlidoul 
2009; Zav ala et alJ 120091 see also llbarra et al.ll2009rfor 

the c ase of decaying dark matter, and Izhang fc SigH 
120081 for dark-matter-induced anisotropics in radio emis- 
sion). When coupled with the energy spectrum, the 
angular power spectrum of the IGRB could be used 
to identify the presence of a da rk matter contribution 
(|Siegal-Gaskins fc Pavlidm][2009t ) . 

In this study, we examine the prospects for using 
modulations in the anisotropy energy spectrum to de- 
tect a signal from Galactic dark matter in the IGRB. 
In §Section [2l we outline the models we adopt for the 
spectral and anisotropy properties of blazars, which are 
taken to be the dominant astrophysical contributor to 
the IGRB, and dark matter. In §Section[3l we determine 
the detectability of a modulation in the anisotropy en- 
ergy spectrum due to a dark matter annihilation signal, 
explore the sensitivity of this technique to variations in 
the blazar and dark matter models, and identify the re- 
gions of the dark matter parameter space accessible to 
Fermi. We discuss the results of this detectability analy- 
sis and its implications for indirect dark matter searches 
with Fermi in §Section|31 

2. THE BLAZAR AND DARK MATTER MODELS 

2.1. IGRB Observables 

We use three metrics to identify a dark matter signal 
in the IGRB: the intensity energy spectrum, the angular 
power spectrum, and the anisotropy energy spectrum. 
The intensity spectrum is the differential photon inten- 
sity / as a function of energy, in units of photons per 
area per time per solid angle per energy. If we expand 
a map of intensity fluctuations SI in units of the mean 
intensity (/), 61 = {I — (/))/ (/), in the basis of spherical 
harmonics, we obtain a set of coefficients aem- The an- 
gular power spectrum Ci = (|a£mp) characterizes these 
fluctuations as a function of angular scale by giving the 
variance of the aim at each multipole i. The anisotropy 
energy spectrum Ci{E) is the angular power spectrum 
evaluated at a specific multipole £ as a function of en- 
ergy E. 

In the case where the total intensity, /tot, is the sum 
of two uncorrelated components (e.g., an extragalactic 
contribution /eg and a Galactic dark matter component 
^dm), the angular power spectrum of the total signal is 
given by 



r2 f^E 



(1) 



where /eg = he /hat and /dm /dm //tot are the frac- 
tional contributions to the total intensity from extra- 
galactic sources and dark matter, respectively. The en- 
ergy dependence of /eg and /dm make the anisotropy 
energy spectrum Ci{E) an energy- dependent function. 

Observation of a feature in the anisotropy energy spec- 
trum of the IGRB would strongly indicate a contribution 
from at least two source classes. On the other hand, if the 
anisotropy energy spectrum is observed to be consistent 
with a constant value at all energies, this would suggest 
that the anisotropy and intensity of the IGRB are domi- 
nated by a single source population at all energies, rather 
than several source populations. An example intensity 



energy spectrum and corresponding anisotropy energy 
spectrum for a scenario in which unresolved blazars and 
Galactic dark matter substructure are the main contrib- 
utors to the IGRB is shown in Figure[Il The dark matter 
contribution is subdominant in the intensity energy spec- 
trum but produces a measurable dip in the anisotropy 
energy spectrum. The error bars represent the l-tr sta- 
tistical uncertainties in the measurement. The calcula- 
tion of the measurement errors and the criteria used to 
determine the energy binning and detectability of a mod- 
ulation are discussed in S Section |3JJ 
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Fig. 1. — Example intensity energy spectrum {top) and 
anisotropy energy spectrum {bottom) for a scenario in which the 
IGRB is composed of emission from Galactic dark matter substruc- 
ture {red dashed lines) and unresolved blazars {blue dot-dashed 
lines) . The total intensity and anisotropy energy spectra are shown 
with solid black lines. The proposed technique tests the measured 
anisotropy energy spectrum for consistency with a constant value. 
See i)Section l3.1l for details. 

We consider a scenario in which the IGRB at ener- 
gies > 1 GeV is composed predominantly of emission 
from unresolved blazars and dark matter annihilation 
in Galactic substructure, as in the example shown in 
Figure [TJ We do not include the contribution of the 
smooth dark matter halo, since it is negligible com- 
pared to the substructure signal at moderate to high 
Galactic latitudes, where the Galactic diffuse emission 
starts to fall off and the I GRB is mos t easily measured 
(e.g., iKuhlen et all 120081: lAndol |2009[) . For simphcity, 
we also neglect the signal from extragalactic dark mat- 
ter, based on predictions that the Galactic substructure 
signal will exceed this em ission (e.g. Hooper fc S erpicol 
l2007t iFtarnasa et aril2009t ). although we note that recent 
work has found that by extrapolating the results of some 
numerical simulations the extragalactic dark matter sig- 
nal can be dominant depending on the structural prop- 
erties assumed for halos and subhalos ([Fieri et al..,2009.) . 
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The Galactic diffuse emission is a major component of 
the gamma-ray sky even at high latitudes and high en- 
ergies, and so we include a treatment of potential uncer- 
tainties from contamination of the IGRB by the Galac- 
tic diffuse emission. Although we do not expect sig- 
nificant anisotropy in this component on 1° angular 
scales (which corresponds to a multipole i ^ 100 as used 
in our analysis), its presence adds statistical noise, and 
residual contamination from this component would re- 
duce the amplitude of the IGRB anisotropy. We have 
approximated the effect of the additional noise from this 
component by taking N\^/Ns = 5, where Ng and A'b are 
the number of signal and background photons, respec- 
tively, observed in the usable sky r egion. We a l so re - 
fer the reader to the recent study of iCuoco et al.l (|20100 
which considers the effect of Galactic foregrounds on the 
sensitivity of Fermi to anisotropics. 

In this study, we choose to focus on a multipole £ — 
100 as this value corresponds to an angular scale that is 
accessible for analysis with Fermi data and also unlikely 
to be heavily contaminated by Galactic diffuse emission. 
We have not evaluated the optimal value of £ for this 
analysis, and expect that different choices of the £ at 
which the anisotropy energy spectrum is evaluated could 
increase or decrease the sensitivity of our method slightly. 

2.2. Spectral and anisotropy properties of blazars 

The intensity spectrum of the collective unresolved 
blazar emission at high energies is still poorly con- 
strained both in shape and amplitude. We model the 
in tensity energy spec t rum of th e blazar emiss i on as 
in lPavlidou fc Venters! pOOl and iPavlidou et al.l (|2008l 
their Appendix C) with the addition of a simple recipe for 
treating gamma-ray attenuation due to the extragalac- 
tic background light (EEL) at high energies. The shape 
of the observed cumulative blazar spectrum depends on 
the blazar spectral index distribution, which in this for- 
malism is parameterized by the mean spectral index of 
blazars ao and the spread in spectral index ag- Given 
a set of observed blazars, a likelihood analysis can be 
performed to find the ma ximum likelihood ao and ctq as 
detailed in IVenters fc Pavlidou (2007 ). Here, we derive 
these parameters from the sample of fiat-spectrum ra- 
dio quasars (FSRQs) in the J'ermi- Large Area Telescop e 
(Fermi-LAT) bright source listH (lAbdo et al.l I2009bn . 
since FSRQs are likely to dominate over BL Lac objects 
in their contribution to the IGRB (.Derm er 2007). We pa- 
rameterize the normalization of the collective blazar in- 
tensity spectrum by its value at 0.26 GeV, Iq = /eg(0.26 
GeV). 

Above a few tens of GeV, extragalactic gamma rays 
start suffering attenuation due to interactions with the 
EBL; for photons in the LAT energy range the relevant 
EBL photons are primarily those of ultraviolet and op- 
tical wavelengths. EBL attenuation induces a suppres- 
sion feature in the cumulative blazar spectrum, which 
can be calculated by integrating the luminosity function, 
properly accounting for the gamma-ray photon optical 

® We caution that faint gamma-ray sources with hard spectra 
are easier to detect than softer sources, and this might bias the 
derived distribution of spectral indices toward harder values if the 
sample used is not flux limited, which is indeed a concern with the 
LAT bright source list. However, especially for the brighter FSRQs 
that we use here, the effect is expected to be small. 



depth given a specific EBL model. However, there are 
large uncertainties both in EBL models and the blazar 
luminosi ty function which p ropagate to the attenuation 
feature (j Venters et al.l [20091 ) . so this calculation would 
unnecessarily complicate our blazar IGRB component 
model in the context of our simple analysis, without 
adding any robust information. Instead, we treat EBL 
attenuation in an approximate manner: we assume that 
most blazar emission is generated by FSRQs at a red- 
shift of zq = 1.0 (a reasonable approximation given the 
curr ently best-available luminosity function models; see, 
e.g. . IVenters et al.]|20 09'). and we accou nt for attenuation 
using the approximate expression in iHoriuchi fc Andol 
(|2006D . We neglect the effects of cascading since we ex- 
pect this to be a minor correction for the models we 
consider here, which have steep mean spectra (of slope 
^ 2.4), so the amount of energy available for cascade 
emiss ion is never do minant over the energy in direct emis- 
sion (IVentersll2009l) . 

Our reference model for the intensity spectrum of the 
blazar IGRB component is shown in Figure [2] For com- 
parison, we overplot a power law / oc consistent 
with Fermi m easurements of the IGRB between 200 MeV 
and 50 GeV (|Abdo et al.ll2010cf) . The reference blazar 
model corresponds to ao = 2.39 and ao = 0.14, the 
maximum likelihood values of the spectral parameters, 
with a normalization at 0.26 GeV of /q = 1.5 x 10~^ 
GeV~^ cm~^ sr~^ to compose the majority of the 
diffuse background between 0.5 and 1 GeV. The depen- 
dence of the detectability of a dark matter feature in the 
anisotropy energy spectrum on the spectral parameters 
and intensity normalization of our blazar model is ex- 
plored in ijSection 13.31 throughout the rest of this study 
we adopt the reference blazar model defined here. 

Here we focus on the anisotropy energy spectrum at 
£ = 100, and therefore C^qq for unresolved blazars is an 
input to our calculation. Unfortunately, the exact ampli- 
tude of the angular power spectrum of the blazar compo- 
nent is dependent on the blazar luminosity function, and 
it is therefore nontrivial to obtain a self-consistent set of 
lo and blazar C^j^qq. The a ngular powe r spect rum of 
blazars has been calculated b ylAndo et al.1 (I2007D. based 
on the lun iinosit y functions of lNarumoto &: Totanil (|2006l 
hereafter ■ lNT06f ) . However, the slope of the fiux distribu- 
tion (logiV-logS') at low fiuxes for Fermi blazars appears 
to be shallower than the NT06 predictions, implying a 
higher anisotropy than the Ando et al. (2007) prediction 
for Fermi. In addition, the level of the Fermi IGRB has 
not been appreciably reduced (at least for energies < 1 
GeV) compared to the determination by Strong et al. 
(2004) for EGRET data. We therefore use the (higher) 
value of C^^QQ calculated by Ando et al. (2007) for the 
EGRET diffuse emission as our reference value. How- 
ever, we stress that an accurate, self-consistent deter- 
mination of the blazar cumulative unresolved emission 
amplitude and anisotropy requires a detailed calculation 
based on luminosity functions derived from Fermi data. 
While the techniques discussed here are robust regard- 
less of the value of C^^qq, the sensitivity curves derived 
could shift up or down depending on how this value dif- 
fers from our assumptions. 

2.3. Spectral and anisotropy properties of dark matter 
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Properties of the dark matter particle x such as its 
mass m^, pair-annihilation cross-section {av), and domi- 
nant annihilation channels, in conjunction with the prop- 
erties of the dark matter subhalo population, will deter- 
mine the shape, amplitude, and location of a dark matter 
feature in the anisotropy energy spectrum of the IGRB. 
The differential photon intensity / in a direction tp from 
dark matter annihilation is 



1 dN (av) 
4^ dE 2m2 



(2) 



where dN/dE is the photon spectrum per annihilation, 
and p{s) is the dark matter density in the line-of-sight 
direction ip at a distance s. The "astrophysical" factor, 
which encodes the mass distribution of the dark matter in 
subhalos and is given above by the line-of-sight integral of 
the dark matter density squared, is set by the structural 
properties of the subhalo population. 

We consider the continuum photon spectra for two 
benchmark annihilation channels: XX ~^ which pro- 
duces a relatively soft spectrum, and xX ~^ t^t~ , which 
produces a harder spectrum. We parameterize the inten- 
sity energy spectr a for these channels us ing the analytic 
approximations of lFornengo et all ()2004[) and show them 
for comparison in Figure [3] Although we defer studying 
the case of gamma-ray lines (xX ~^ 77 or xX ~^ ^'^l) 
to future work, we note that the anisotropy energy spec- 
trum may be highly sensitive to these annihilation chan- 
nels since the line feature in the intensity spectrum could 
lead to a dramatic variation in fumiE) against the col- 
lective emission s pectr um of astrophysical sources (see 
iZhang fc Beacoinl f200l. 

We adopt the fiducial model Al of lAndol (|2009f ) to 
self-consistently determine the amplitude of the angular 
power spectrum and the astrophysical factor entering in 
the normalization of the dark matter intensity. Model Al 
extrapolates the results of recent numerical simulations 
to a minimum subhalo mass M^-^ = 10~^ Mq, using a 
mass function slope a = 1.9, which leads to a fraction 
/ = 0.2 of the total halo mass in substructure for this 
choice of M^i^ and normalization of the mass function. 

3. DETECTABILITY OF A DARK MATTER 
ANNIHILATION FEATURE 

Having defined the models for blazars and dark matter 
which will comprise the two-component scenarios consid- 
ered in this study, we now determine in what cases the 
presence of the dark matter component can produce a de- 
tectable modulation in the anisotropy energy spectrum. 

3.1. Detectability Criterion 

In order to quantitatively assess the sensitivity of the 
anisotropy energy spectrum as a tool for indirect dark 
matter detection, we need to define a uniform detectabil- 
ity criterion. We do this by testing the hypothesis that, 
in each scenario we examine, the anisotropy energy spec- 
trum at energies > 0.5 GeV is consistent with a constant 
value, equal to the weighted average fi of all energy bins 
in which the measured Ce is at least 2-a from zero. If this 
hypothesis is rejected by a x^ test at a 95% confidence 
level (CL), then we consider the particular dark matter 
model detectable. 

To evaluate the sensitivity of this technique using 
Fermi data, we assume that a fraction of the sky /sky 
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Fig. 2. — Reference model for the blazar contribution to the 
IGRB (black dot-dashed line), with parameters oq = 2.39, (tq = 
0.14, /q = 1-5 X 10~5 GeV~i cm'^ g-i sr-\ and zo = 1. For 
comparison we show a power law I oc (red solid line), con- 

sistent with Fermi measuremen ts of the IGRB between 200 MeV 
and 50 GeV IIAbdo et al.ll2010d) . 
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Fig. 3. — Continuum photon spectra for a m-^ = 500 GeV dark 
matter particle annihilating to bb {purple solid line) and to t~^t~ 
(blue dashed line). 



of 0.75 will yield data usable for this analysis, since we 
expect that contamination by Galactic foregrounds will 
be substantial in the inner Galaxy. We estimate the num- 
ber of signal (dark matter plus blazar) photons that 
will be detected in the usable fraction of the sky in one 
or five years of all-sky observation time. Fermi operates 
primarily in sky-scanning mode, resulting in fairly uni- 
form exposure across the sky; hence, we assume uniform 
exposure. We take the field of view = 2.4 sr, and ap- 
proximate the effective area Ace as a function of energy 
according to the reported Fermi-hAT performanceQ. 

Bins in energy are logarithmic and correspond to an 
increase in energy of a factor of 2, with the additional 
constraint that at least 2 x 10** signal (dark matter plus 
blazar) photons are detected in each bin. This is esti- 
mated by convolving the total intensity spectra of the 
dark matter and blazar model being tested with the ef- 
fective area of Fermi then multiplying by the observation 
time. The 1-a statistical uncertainties in each energy bin 

^ |http: //www-glast ■ slac ■ Stanford. edu/software/IS/glast_lat_performaiic 
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are given by (see, e.g.. lKnox|[l995l ) 



AC? = 



(2^+l)A£/,ky 



(3) 



where C| is the angular power spectrum of the signal (in 
our case, C|°') and Wg — exp(— i?^cr^/2) is the window 
function of a Gaussian beam of width db, where db is 
the energy-dependent LAT angular resolution, which we 
approximate by a parameterization of the reported LAT 
performance. We use multipole bins of A£ — 10. The 
noise power spectrum Cn = (47r/sky/iVs)(l + {N^/Ns)) is 
the sum of the Poisson noise of the signal (here blazars 
and dark matter) and the background (e.g., diffuse emis- 
sion from Galactic cosmic ray interactions and instru- 
ment noise). The ratio of the number of background 
photons to the n umber of signal pho tons N\,/Ns is taken 
to be 5 (see, e.g.. lAbdo etal1l2009c| ). 

For each scenario, defined by a specific blazar and dark 
matter model, we test the consistency of the observable 
anisotropy energy spectrum with a constant value by ap- 
plying a test. The standard statistic is: 



E 



(4) 



where O is the measured value of Ci at each energy bin 
and a = AC| is the 1-a error on each poinfl We have 
attempted to appropriately estimate the measurement 
uncertainties taking into account the expected signals, 
backgrounds, and instrumental capabilities; however, we 
caution that underestimating the measurement uncer- 
tainties would result in an overestimate of the sensitivity 
of the technique and could lead to spurious detections. 
To apply this technique to LAT data, a careful determi- 
nation of the measurement uncertainties will be essential. 

In summary, we first compute the total intensity energy 
spectrum of the blazar model combined with the dark 
matter model being tested as would be seen by Fermi 
by using the spectra defined in §Section[2]and the Fermi 
observational parameters described in this section. Us- 
ing this intensity energy spectrum and the Ci values for 
blazars and dark matter discussed in §Section[21 we com- 
pute the anisotropy energy spectrum via Equation ([1} to 
determine the Ci value which would be observed in each 
energy bin, with errors given by Equation ([3|). Using 
Equation we perform a test on this spectrum to 
determine if the anisotropy energy spectrum is consistent 
with a constant value at the 95% CL. If it is not consis- 
tent with a constant, then we declare the modulation in 
the anisotropy energy spectrum due to the dark matter 
component to be detectable. 

3.2. Sensitivity to dark matter models 

We now examine the sensitivity of this technique to 
dark matter particle properties. In FigureHl we show the 

® The uncertainties used in calculating the weighted average of 
all bins arc calculated using the value of C( measured in each bin in 
Equation JSj. However, the uncertainties that enter Equation (|4]l 
are calculated by using the weighted average value of Ci in Equa- 
tion ^ for all bins, as this is now the hypothesis that is being 
tested. This choice, in combination with the fact that for £ = 100 
the probability distribution of Ci asymptotically approaches a nor- 
mal distribution, renders the uncertainties approximately Gaus- 
sian, and the CLs of the standard analysis apply. 



detectable region of the m-^-{av) parameter space for one 
and five years of Fermi observations assuming our refer- 
ence blazar model and annihilation to bb or t~^t~^ . Mod- 
els above the lines for each annihilation channel are de- 
tectable for the test criterion defined in Section [XT] Us- 
ing this criterion, we find that annihilation cross-sections 
of order 3 x 10~^^ cm^ s~^, the expected value for a ther- 
mal relic, are within the reach of Fermi for the t'^t~ 
channel for masses up to ^ 100 GeV (50 GeV) with five 
years (one year) of observations, and cross sections within 
a factor of 10 of thermal are detectable for a large range 
of masses (between ^ 10 and ~ 400 GeV) within five 
years for this annihilation channel. For annihilation to 
bb, cross sections within a factor of 10 of thermal are 
detectable in five years for the entire range of masses 
considered here, 10 GeV < < 1 TeV. It is notable 
that this technique is sensitive to a viable region of the 
parameter space; existing constraints robustly exclude 
total annihilation cross sections larger tha n ^ 10"^'^ cm'^ 
s~^ f or the mass range considered here (jBeacom et al.l 
l2007f ). Reionization constraints from the observed opti- 
cal depth of the universe more strongly limit the cross- 
section for these annihilation channels to (erf) < 10~^* 
100 GeV (e.g., fSlatver et al.l 12001 



< 



IHtitsi et al.l[200aiCirelh"etaII[2009l) . 

The subtle bump-like features in Figure SI present for 
both annihilation channels, are a result of the specific cri- 
terion we chose for the statistical test, and not an indica- 
tion of a rapid change in the properties of the dark mat- 
ter model at those masses. Specifically, the energy range 
and multipole we selected for this analysis lead to weaker 
sensitivity in these regions. This can be understood in 



the bb one- year case by noting that for 100 



< 



< 200 



GeV, this channel produces a large number of photons at 
energies of a few GeV, where Fermi can generally make a 
precise measurement of the anisotropy. As the blazar and 
dark matter intensity spectra fall above this energy, the 
number of energy bins in which a precise measurement 
of the angular power spectrum can be made is limited 
due to low photon statistics. In this situation, a devia- 
tion in the anisotropy energy spectrum (i.e., the transi- 
tion at high energies to the anisotropy being dominated 
by blazars) is difficult to detect by our criterion using a 
small number of points which may have relatively large 
measurement uncertainties. The features in the t'^t~ 
sensitivity curves are present for similar reasons. Alter- 
native parameters for the statistical test (e.g., a different 
multipole or energy range) would modify these features. 
The sensitivity of this technique to these models may 
improve by optimizing the detection criterion of the test 
considered here, or by the application of more tailored 
likelihood analysis rn ethods to test these scenarios (e.g., 
iDodelson et al.|[2009D : however, such a detailed study is 
beyond the scope of this work. 

As an indication of the relative level of dark matter 
contribution to the IGRB in our models, we show in Fig- 
ure [5] the annihilation cross section (av) at which the 
dark matter intensity exceeds that of the reference blazar 
model at any energy E > 0.5 GeV as a function of dark 
matter mass for the two annihilation channels con- 
sidered {solid curves). As expected, for large dark matter 
masses (m^ > 300 GeV), the intensity of the reference 
blazar model is exceeded even at fairly small annihila- 
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Fig. 4. — Sensitivity of Fermi to modulations in the anisotropy 
energy spectrum due to dark matter ann ihila tion in substructure 
for the test criterion specified in Sectioji 13.11 The minimum de- 
tectable annihilation cross section to bb (purple) or r+r" [light 
blue) is shown assuming one year (solid) or five years (dashed) of 
observations. Dark matter models above the curves produce a fea- 
ture detectable at 95% CL. 
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Fig. 5. — Annihilation cross-section (av) as a function of dark 
matter mass m-^ above which the dark matter intensity exceeds the 
intensity of the reference blazar model at any energy above 0.5 GeV 
(solid curves), for annihilation to bb (purple) and t^t~ (light blue). 
Sensitivity curves from Figure|4]for five years of observation assum- 
ing the reference blazar model are shown for comparison (dashed 
curves). The dark matter signal is subdominant in the intensity 
energy spectrum at energies above 0.5 GeV for models below the 
solid curves. 

tion cross sections for the t+t^ channel, due to the ex- 
ponential cut-ofF of the blazar intensity above a few tens 
of GeV and the relatively hard_continuuin spectrum of 
this channel. Annihilation to bb produces a softer con- 
tinuum spectrum, and consequently larger annihilation 
cross sections are required to exceed the blazar intensity 
for the mass range considered here than in the t~^t~ case. 
The detectable region of the parameter space from Fig- 
ure|3]for five years of observation is shown for comparison 
{dashed curves). Comparing the regions indicated by the 
solid and dashed curves, it is clear that the dark matter 
component does not need to be dominant in the intensity 
to produce a signal detectable by this technique. 

In Figure |6l the ratio of the dark matter to blazar 
intensity at 1, 10, and 100 GeV is shown for dark mat- 
ter models detectable within five years at 95% CL. The 
normalization of the dark matter intensity at each mass 
is set by the annihilation cross section corresponding to 




10 100 1000 
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Fig. 6. — Ratio of dark matter to blazar intensity /dm/^eg at 
three reference energies as a function of dark matter mass for 
detectable models (annihilation cross sections) corresponding to 
the five-year sensitivity curves shown in Figure [J] Relative inten- 
sity of the dark matter signal is shown at i? = 1, 10, and 100 GeV 
(as labeled) for annihilation to bb (purple) and r+r~ (light blue). 
The E = 100 GeV curves do not extend to masses below 100 GeV 
because no dark matter annihilation photons arc produced at 100 
GeV for < 100 GeV. 

the five-year sensitivity curves shown in Figure 3] for each 
annihilation channel. A dark matter intensity < 10% of 
the blazar intensity a.t E — 1 and 10 GeV is sulRcient 
to produce a detectable feature in the anisotropy energy 
spectrum for the mass range considered here (10 GeV - 
1 TeV). At E = 100 GeV, detectable models contribute 
a larger intensity relative to the blazar intensity, up to 
~ 40% at = 1 TeV for annihilation to bb and up to 
~ 200% at = 1 TeV for annihilation to r+r^. How- 
ever, we note that our reference blazar model does not 
account for most of the Fermi measurement of the IGRB 
at E = 100 GeV (see Figure [21 assuming the power- law 
IGRB can be extended to i? = 100 GeV), so dark mat- 
ter models that exceed our blazar model intensity at this 
energy by a factor of ~ 2 are still acceptable given the 
measured IGRB level. 

3.3. Dependence on blazar model parameters 

We now examine how the detectability of a dark mat- 
ter feature in the anisotropy energy spectrum depends on 
the overall normalization and the spectral shape of the 
blazar component of the IGRB. We show in Figure [7] the 
sensitivity of this technique to variations within the l-cr 
uncertainty in our determination of the blazar spectral 
parameters ao and ctq (^shaded bands) and variations in 
the blazar normalization, assuming five years of obser- 
vations. The sensitivity for the blazar intensity normal- 
ized to 92% (dotted lines) and 30% {dashed lines) of the 
reference model normalization, assuming the maximum 
likelihood spectral parameters, is also shown. These re- 
sults indicate that the uncertainty in the test sensitivity 
introduced by uncertainties in ao and cto is small com- 
pared with the range of detectable signals. Reducing the 
blazar normalization to the extent shown here improves 
our test sensitivity slightly, although it is important to 
keep in mind that our simple 2-component model would 
become unrealistic for very small blazar normalizations, 
since if blazars are a subdominant component of the mea- 
sured IGRB we expect other source classes, not consid- 
ered here, to contribute significantly. Other scenarios 
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for astrophysical sources of the IGRB will be addressed 
in a future publication. We note, however, that in the 
likely scenario where a population of numerous but in- 
dividually faint sources (such as star-forming galaxies or 
quiescent blazars) contributes a substantial fraction of 
the IGRB, the angular power at ^ '-^ 100 could still be 
dominated by bright blazars, since such an additional 
population would ap pear very isotropic in co mparison 
on small scales (e.g., lAndo fc Pavlidoul |2009() . In this 
case, the additional component would behave as an al- 
most isotropic "background" on top of our bright blazar 
plus dark matter "signal" . 




10 100 1000 

[GeV] 



Fig. 7. — Minimum value of the annihilation cross section for 
which a dark matter annihilation signal is detectable at 95% CL in 
five years of observation with Fermi as a function of , for anni- 
hilation to bb {purple) or t~^t~ {blue). The height of the band at 
each dark matter mass represents the uncertainty in detectability 
entering through the l-cr uncertainties in the blazar spectral pa- 
rameters ao and erg, assuming the reference blazar intensity nor- 
malization; the solid line within each band shows the sensitivity 
obtained for the maximum likelihood values of the spectral param- 
eters. The sensitivity curves for the blazar intensity normalized to 
92% {dotted lines) and 30% {dashed lines) of the reference model 
normalization, assuming the maximum likelihood spectral param- 
eters, are shown for comparison. 

4. DISCUSSION 

We have examined the potential of the anisotropy en- 
ergy spectrum of the IGRB to identify a dark matter an- 
nihilation signal for the case that the main contributors 
to the IGRB are unresolved blazars and Galactic dark 
matter subhalos. We have evaluated the detectability of 
a dark matter signature as a function of the properties of 
the dark matter particle and of the blazar population us- 
ing a model-independent detectability criterion, namely 
the consistency of the anisotropy energy spectrum with 
a constant value over varying energies. 

We find that modulations in the anisotropy energy 
spectrum are a sensitive probe of dark matter models. 
Using our reference blazar model and assuming annihi- 
lation to T+r~, Fermi could detect in one year a fea- 
ture from dark matter for annihilation cross sections as 
small as that expected for a thermal relic for particle 
masses below ^ 50 GeV. In five years, cross sections 
within a factor of 10 of thermal produce a detectable 
feature for 10 < < 1 TeV for annihilation to bb, and 
10 ^ ™x ^ 400 GeV for annihilation to r+r". The 
dependence of the test sensitivity on the blazar model 



parameters (spectral shape and overall normalization) is 
minimal and does not affect our conclusions qualitatively. 

We have treated the simple case where only blazars 
and Galactic dark matter have a significant contribution 
to the IGRB. In reality, other classes of sources, such 
as extragalactic dark matter, normal star-forming and 
starburst galaxies, and galaxy clusters may also be sig- 
nificant IGRB components. We plan to study the effect 
of these additional contributions on Galactic dark matter 
detectability in a future publication. 

Recent results from Fermi (Abd o et al.l l2010d( ) in- 
dicate that the total contribution of sources of the 
same classes as those comprising the majority of Fermi- 
resolved sources (such as bright blazars) may comprise 
as little as 10%-30% of the total IGRB intensity. In 
that case, the additional IGRB photons must originate 
in truly diffuse processes or unresolved objects that are 
underrepresented in the set of already-resolved Fermi 
sourc es. Examples of the latt er are quiescent blazars 
e.g.. iStecker fc SalamonI |1996[) . star- forming galaxies 
Fields et al. I i2oion . and dark matter annihilation in 
Galactic or extragalactic structures. However, as dis- 
cussed in § Section 3.3, it is likely that even if the in- 
tensity is not dominated by bright blazars, the angular 
power at small scales is: numerous but individually faint 
sources will contribute very little to the total anisotropy. 

There is at least one guaranteed contribution to the 
IGRB which is expected to be similar in spectral shape 
with dark matter annihilating to bb (although peak- 
ing at lower energies than most plausible dark mat- 
ter candidates): millisecond pulsars, which constitute 
a kn own Fermi source population with resolved mem- 
bers (jAbdo et al.ll2009al ). and which have a spatial dis- 
tribution extending to high Galactic latitudes (see, e.g., 
Faucher-Giguere & Locb 2010). However, their contri- 
bution would be limited to the lower end of the energy 
range we have considered, as their spectra generally cut 
off exponentially above ~ few GeV. In addition, since 
the MSPs are correlated with the Galactic plane, the 
large-scale morphology of their contribution to the dif- 
fuse emission would show a gradient away from the plane, 
which could be used to identify its origin. Finally, the 
spectral signature of millisecond pulsars is fairly well con- 
strained, so in a more detailed study the modulation of 
the anisotropy energy spectrum due to the presence of 
millisecond pulsars can be explicitly modeled and taken 
into account. 

The sensitivity estimated in this work is based on a 
generic test for the presence of any detectable deviation 
from energy invariance in the anisotropy energy spec- 
trum; the results obtained depend on the parameters 
defining the statistical test, such as the choice of mul- 
tipolc and energy bins. To straightforwardly define our 
statistical test, we have not fully optimized the param- 
eters of the test for each set of dark matter parameters. 
This is reflected in the features in the sensitivity curves 
shown in Figure [H which arise from the parameters cho- 
sen for our statistical test, and can change qualitatively 
for different test parameters. We expect that instead 
testing the likelihood of a specific dark matter and blazar 
model producing observed data will extend the sensitiv- 
ity of this technique, and defer that analysis to future 
work. 

The anisotropy and overall intensity in the dark mat- 
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ter subhalo gamma-ray signal was calculated in a manner 
as self-consist ent as possib le based on one of the fidu- 
cial models of lAndd (|2009[ ). However, this model does 
not take into accou nt the radial depend ence of subhalo 
concentrations fe.g-. lKuhlen et alJ l2008') which may re- 
sult in increased anisotropy and intensity from subha- 
los. In addition, a steeper subhalo mass-luminosity re- 
lation than employed in this model, especially for sub- 
halos be low ^ 10^ M^, is p redicted from the simula- 
tions of iKuhlen et alJ (|2008l ). Using this steeper rela- 
tion would decrease the contribution from low-mass sub- 
halos to the intensity, which would likely increase the 
overall amplitude of the angular power spectrum from 
Galactic substructure. The anisotropy of the blazar sig- 
nal wa s calculated based on the luminosity functions of 
NT06I. However, first results from Fermi (jAbdo et alJ 
2010d( ) may indicate that the nu mber de nsity of blazars 
is lower than the predictions of INT061 which, if con- 
firmed, would imply that the anisotropy in the blazar 
signal may also be higher than the values considered 
here. These effects would likely lead to improved de- 
tectability and increased constraining power since we are 
generally in a low-statistics regime, and thus larger am- 



plitude anisotropics can be measured more easily. 
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